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We present preliminary results on two new two-particle correlation analyses of Au+Au collisions 
at -y/ijviv" = 13QGeV performed by the STAR collaboration. Two-pion interferometry with respect 
to the reaction plane shows that the pion source can be described as an ellipse extended out of the 
reaction plane. Analysing the pion-kaon correlation functions, we show that, on average, pions and 
kaons are not emitted at the same position and/or time. Both measurements are interpreted in the 
so-called blast wave model framework whose main feature is a strong flow. 



I. INTRODUCTION 



Two-particle interferometry is a sensitive probe of the 
space-time geometry of the particle emitting source. Us- 
ing this technique, pion source sizes have been measured 
at every available relativistic heavy ion beam energy 
which led to the excitation function reported in |IJ . Two- 
particle correlation also offers a large variety of observ- 
ables combining different particle species. For example, 
two-kaon and two-proton correlation functions have been 
measured at the CERN SPS j|, [|. I n this paper, we 
present preliminary results from two-particle correlation 
analyses of Au+Au collisions at ^/snn = 130 GeV pro- 
duced by the Relativistic Heavy Ion Collider (RHIC) at 
Brookhaven National Laboratory. 

Two-particle correlations are studied by construct- 
ing correlation functions : C 2 (k*) — A(k*)/B(k*) with 

k* = Pl ~ P2 the momentum of one particle in the rest 
frame of pair, A(k*) the distribution of k* for pairs 
of particles from the same event, and B(k*) the k* 
distribution for pairs of particle from different events. 
Note that for pairs of particles ha ving the same mass, 
Qinv = \J (pi - P2) 2 - {Ei - E 2 ) 2 = 2k*, where Q inv 
is the invariant relative momentum. For pairs of par- 
ticles having different masses, this relationship doesn't 
hold and k* is the only relevant variable. 

Transverse mass spectra Q] and elliptic flow analysis || 
suggest that there is a strong collective motion of the par- 
ticles in the Au+Au collisions at ^snn = 130 GeV. Such 
flow introduces a strong correlation between the position 
and the momentum of the emitted particles, which re- 
duces the apparent source size. It could explain why the 
analysis of two-pion correlation functions published by 
STAR jjj shows no anomalously large pion source size 
which would have indicated the formation of the quark 
gluon plasma. Transverse flow must be accounted for 
in order to interpret the source size extracted from two- 
particle correlation analysis. We will introduce the so 
called "extended blast wave model" in order to investi- 
gate the effect of flow. 

We present two new analyses that are potentially more 
sensitive to flow effects. They both provide new indepen- 
dent insights on how the system looks at kinetic freeze- 
out, i.e. when the particles decouple. 



II. THE EXTENDED BLAST WAVE MODEL 

The extended blast wave model allows one to combine 
transverse mass spectra, elliptic flow and two-particle 
correlation analysis within a single framework. It is a 
simple parametrization of the system at kinetic freeze- 
out in terms of temperature, transverse flow, and source 
transverse radius. This approach focuses on the trans- 
verse plane, ignoring the longitudinal direction, i.e., the 
direction parallel to the beam. The blast wave model was 
introduced in to interpret the identified particle el- 
liptic flow from STAR. In this model, particles are emit- 
ted from a infinitely thin shell. To reproduce the pion 
source size measured, the model had to be extended to a 
filled cylinder. Such a parametrisation had been used for 
several years to interpret transverse mass spectra || |?j] 
assuming a anzymuthally symetric system. 

In this extended blast wave framework, the probability 
of emitting a particle at a given space-time X and energy- 
momentum P is : 

f( x , p ) = f{r,0s,t,PT,<f> P ,m) 
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This equation can be understood separating it in 3 parts: 

- The step function, 8 confines the system inside a filled 
ellipse. 4> s ((f>p) is the spatial (momentum) azymuthal 
angle ; (f> s = for particle emitted in the reaction plane. 
R x and R y are the maximum radii of the system in plane 
and out of plane, respsectively. To interpret the identified 
particle elliptic flow Q , the variable s 2 was introduced as 
a parameter that quantifies the spatial anisotropy of the 
particle emission in non-central collisions. We find, when 

1 3 — 1 

describing the system by an ellipse, that s 2 ~ ^ with 
r\ = R y /R x . When s 2 > 0, Ry > R x which means that 
the source is described by an ellipse extended out of the 
reaction plane. 

- The second term of the equation express the hydro-like 
behavior of the system. The expression was derived in 
|Q where the following functions are used : a(r, 4> s ,pt) = 
^ smh(p(r, </>„)) and f3(r, cj> s ,m T ) = ^ cosh(p(r, <j> a )). p 
is the transverse rapidity boost that the particles feel 
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at a given emission point. In order to describe non- 
azymuthally symetric systems, we have modified its orig- 
inal expression, to include its dependence on the azy- 
muthal angle : p{r,4> s ) = |^-(po + Pa cos(20 s )), where 

f = yj (r ■ sin(4>s)) 2 + r] 2 (r ■ cos(4> s )) 2 is the "elliptical 
radius" . I.e. all points on an ellipse have a constant 
f. p a > means that there is a stronger boost in 
the reaction plane than out of it. Since in an ellipse 
the radial vector and the vector normal to the surface 
of the ellipse are not aligned we introduce the variable 
<j>b — tan~ 1 (taxi(4> s )/ri 2 ). Indeed, in the Cooper Frye 
prescription || that the blast wave model follows, the 
transverse boost that the particles acquire has to be nor- 
mal to the freeze-out surface. 

- The third term of the equation accounts for the dura- 
tion of particle emission using a single parameter r. This 
term is particulary important when interpreting the pion 
source size. It has no effect, however, on the transverse 
mass spectra and elliptic flow. 

In the azymuthally symetric case, this formula simpli- 
fies to : 

f( x , p ) = f(r,0s,t,PT,(l>p,m) 
= G(R-r) ■ 

K x {p{r,m T )) • e «(nPT)-co S (^ s -0 p ) . 
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with p(r) — |po-^- The system is now confined to a 
cylinder of radius R. There are four parameters: the tem- 
perature T, the magnitude of the flow transverse boost 
po, the radius of the cylinder R, and the emission dura- 
tion t. The parameters that quantify the anisotropy of 
the system (s 2 and p a ) vanish. 

From the phase space density, f(X,P), it is straight- 
forward to calculate transverse mass spectra and source 
radii. Both the pion source radii and the pion, kaon and 
proton transverse mass spectra are well reproduced with 
po = 0.6, T = 110 MeV, R = 13 fm and r = 1.5 fm/c. It 
is interesting to notice that to keep the ratio R ou t / Rside 
as a function of px comparable with the values measured 
by STAR[jl| which are close or even below one, the emis- 
sion duration must be on the order of 1 or 2 fm. Indeed, 
since in the usual Bertsch-Pratt parameterization, R ou t 
is the radius component parallel to the transverse mo- 
mentum of the pair (hp) while R S ide is normal to kx, the 
emission duration increases R ou t while leaving R S id e un- 
changed. Such a short emission duration is not currently 
achieved by any realistic microscopic or hydrodynamic 
model. 

The parameters that best fit the identified particle el- 
liptic flow measured by STAR § are p = 0.61 ± 0.05, 
T = 101±24MeV, p a = 0.04±0.01 and s 2 = 0.04±0.02. 
The flow rapidity and temperature are consistent with 
the values obtained from transverse mass spectra and 
pion source size. The system is best described when both 
p a and s 2 differ from zero, i.e. when the system exhibit 
an asymmetry in momentum and space. 



III. PION SOURCE GEOMETRY WITH 
RESPECT TO THE REACTION PLANE 

The identified particle elliptic flow favours a system 
that is asymmetric both in space and momentum. The 
pion source should reflect such an asymmetry. It is thus 
interesting to study the pion source geometry as a func- 
tion of the angle between the particle momenta and the 
reaction plane. 

Such analysis was performed by the E895 collabora- 
tion |j| at the AGS. The important difference between 
STAR and E895 is that directed flow (vi) cannot yet 
be measured by STAR, which means that only the radii 
within the transverse plane are relevant. The available 
radii are R ou t, Rside and the cross term R ou tside- For the 
theoretical aspect of the analysis see ]Tc| . 

Events were selected using a minimum bias trigger 
which covers 85% of the whole centrality range. It is 
mandatory to include peripheral events where the ellip- 
tic flow is maximal. Pion tracks are reconstructed and 
identified in the STAR TPC Jll| in the momentum range: 
0.1 < p T < 0.6 GeV/c and -0.5 < Y < 0.5. The reaction 
plane is reconstructed as described in S . Pion pairs are 
sorted into four different bins depending on the angle be- 
tween the pair momentum and the reaction plane. R ou t, 
Rside and Routside are extracted fitting the four different 
correlation functions constructed with each pair sample. 
Figure | shows the squared radii R 2 out , R 2 ^ de and R 2 outslde 
as a function of the mean angle of the pairs with respect 
to the reaction plane. 

A clear oscillation of the radii is observed. The line in 
the figure is a blast wave model calculation using the pa- 
rameters po = 0.6, T = 100 MeV, p a = 0.05, s 2 = 0.05, 
R = 10 fm, and r = 2 fm/c. These parameters were cho- 
sen so that they are consistent with the ones extracted 
from identified particle elliptic flow and two-pion intcr- 
fcromctry. A non zero value of the space asymmetry 
parameter s 2 is necessary to reproduce the data. When 
s 2 is equal to zero the momentum asymmetry expressed 
by p a is not sufficient to reproduce the oscillation of the 
radii. Thus, since s 2 > 0, this analysis indicates that the 
pion source geometry is an ellipse extended out of the 
reaction plane. 



IV. PION-KAON CORRELATION FUNCTION 

The correlation of non-identical particles can be used 
to study the relative mean separation between the emis- 
sion time and/or position of different particle species |l2] , 
P|. This technique is based on the comparison of two 
different correlation functions. If one assumes that kaons 
and pions are not emitted, on average, at the same time 
or position, for each pion-kaon pair there are two configu- 
rations: the pion and kaon get closer to each other or they 
move away from one another. The correlation between 
the two particles will be stronger in the first configura- 
tion than in the second one. It is shown in [fL2l that 
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FIG. 1: Pion source radii as a function of the angle between 
the pair momenta and the reaction plane. Triangle : data. 
Solid line : blast wave model calculation. 

these two configurations correspond to ~lf pa ir ■ k either 
positive or negative ( ifpair is the velocity of the pair and 
k is the momentum of either the kaon or the pion in the 
rest frame of the pair). Two correlation functions can 
be constructed : C+(| k |) for pairs with v pa il • f? > 
and C_(| 1? |) for 'v'pair • A; < 0. If both correlation 
functions are not identical, it implies that there is a dif- 
ference in the average space-time emission point of pions 
and kaons, or in other words, there is a space-time asym- 
metry in the emission process of pions and kaons. Thus, 
to study whether or not pions and kaons arc emitted at 
the same time we calculate the ratio C_(fc*)/C + (fc*) as 
a function of k* =\ kr \. 

In this analysis, only the most central events are se- 
lected. They represent 12% of the total hadronic cross 
section. Pions and kaons are identified by their specific 
energy loss in the STAR TPC. The acceptance of the pi- 
ons is : 0.1 < p T < 0.6 GeV/c and -0.5 < Y < 0.5, 
and 0.3 < p T < 0.8 GeV/c and -0.5 < Y < 0.5 for 
the kaons. Electrons and positrons are carefully removed 
to avoid contamination of the ir + — K~, and ir~ — K + 
pairs by e + — e~ pairs which are correlated in a non- 
trivial way. Special care is taken to remove two-track 
merging effect since it noticeably influences the 7r + — K + 
and 7T _ — K~ correlation functions. The purity of the 



pion and kaon sample needs to be precisely known in 
order to estimate source size parameters from the cor- 
relation function. Indeed misidentified particles, or sec- 
ondary particles such as pions from K® decay lower the 
correlation strength. To first order, this effect is equiv- 
alent to increasing the source size which also lowers the 
average correlation strength. The purity analysis is still 
under way which prevents us from making any quantita- 
tive statement from this analysis. 

The correlation between the pions and kaons arises 
mainly from their Coulomb attraction if they have op- 
posite charge or repulsion if they have the same charge. 
The 7r + — K + , 7r~ — K~ . 7T + — K~ and ir~ — K + cor- 
relation functions shown in figure |2j exhibit such a be- 
haviour. A direct comparison to the blast wave model 
cannot be performed because these correlation functions 
depend not only on the transverse source size but also on 
the longitudinal one, which is not accounted for in the 
blast wave framework. 
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FIG. 2: 7T+ - K+, IT' - K-, 7T+ - K-, and w~ - K+ 
correlation functions 

The blast wave predictions can however be compared 
to the ratio of correlation functions C_(fc*) over C+(fc*) 
shown in figure |3[ Indeed, the acceptance of each particle 
species is symmetric about mid-rapidity, which cancels 
out any effect that would be due to a difference in the 
emission position along the longitudinal direction. To 
improve the statistics, we have combined the 7r + — K + 
and 7T _ — K~ correlation functions and the 7r + — K~ and 
7r~— K + correlation functions. The ratio C_(fc*)/C + (fc*) 
is significantly different from unity for both like sign and 
unlike sign pairs which leads us to the conclusion that 
pions and kaons are not emitted at the same space or/and 
time. 

The extended blast wave model calculations compare 
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FIG. 3: Ratio of correlation functions C_(fc*)/C+(fc*) 



correlation, which leaves the particles more freedom to 
fill the cylinder volume. Thus, the thermal motion pulls 
the average emission radius away from the edge of the 
cylinder, towards the centre of the system. Kaons, be- 
ing less affected by the thermal motion due to their high 
mass compare to the temperature, are more likely to be 
emitted close to the edge of the cylinder than pions. In 
other words, on average, kaons are emitted further from 
the centre of the source than pions which is why the ra- 
tio C_(fc*)/C+(fc*) calculated in the extended blast wave 
framework, shown on figure ^ deviates from unity. 

The interplay between the temperature and transverse 
flow built in the blast wave model introduces an asym- 
metry in the emission of the pions and kaons that is in 
qualitative agreement with the data. The non-identical 
particle correlation function can be used to constrain the 
temperature and flow parameters of the extended blast 
wave model. However, at this stage, given the statistical 
and systematic uncertainties we only conclude that the 
blast wave model is in agreement with the data when us- 
ing the same parameters that reproduce the transverse 
mass spectra and the pion source radii. 



well with the data. The calculation was done with 
po = 0.6, T = 110 MeV, R = 13 fm, and r = fm/c. 
We showed in section II that, with these parameters, the 
extended blast wave model describes well the transverse 
mass spectra and the pion source size. It is striking to 
notice that the extended blast wave model reproduces 
the data without the need of any extra parameters such 
as one representing the difference between the average 
emission time of pions and kaons. Why does the blast 
wave model predict an asymmetry between the emission 
point of the pions and kaons? To answer this question 
it is important to notice that at low k*, pions and kaons 
have the same velocity but not the same momentum. In- 
deed, at low k* the pion average is equal to 0.1 GeV/c 
while it is equal to 0.4 GeV/c for the kaons. In the blast 
wave model, if the temperature is equal to zero, particles 
are emitted at radii that depend only on their transverse 
velocity not on their mass. In this case, pions and kaons 
would be emitted at the same radii. However, switching 
on the temperature decreases the position-momentum 



V. CONCLUSIONS 

We have shown two new analyses from the STAR ex- 
periment at RHIC. The analysis of the pion source ge- 
ometry with respect to the reaction plane led to the 
conclusion that the source is an ellipse extended out of 
plane. The study of the pion-kaon correlation functions 
allows us to state that < px >= 0.1 GeV/c pions and 
< pt >= 0.4 GeV/c kaons are not emitted at the same 
position or/and time. 

We have interpreted our results in the framework of 
the extended blast wave model. We show that the pion 
source geometry, and the pion-kaon space-time emission 
asymmetry can be described by the same set of parame- 
ters which are also used to interpret the transverse mass 
spectra and the elliptic flow measured by STAR. The 
picture that emerges from this model is a system with a 
strong collective expansion that freezes out in a few fm/c. 
Such a scenario is not currently achieved by any realistic 
models, hydrodynamic or microscopic. 
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